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Abstract. We describe a simple n-dimensional quantum cellular au- 
O^l tomaton (QCA) capable of simulating all others, in that the initial config- 

uration and the forward evolution of any n-dimensional QCA can be en- 
coded within the initial configuration of the intrinsically universal QCA. 
pH Several steps of the intrinsically universal QCA then correspond to one 

O | step of the simulated QCA. The simulation preserves the topology in 

I the sense that each cell of the simulated QCA is encoded as a group of 

adjacent cells in the universal QCA. 

, 1 Introduction 

Cellular automata (CA), first introduced by Von Neumann [3B], consist of an 
array of identical cells, each of which may take one of a finite number of possible 
states. The whole array evolves in discrete time steps by iterating a function 
G. This global evolution G is shift-invariant (it acts everywhere the same) and 
7-H local (information cannot be transmitted faster than some fixed number of cells 

per time step). Because this is a physics-like model of computation [15] . Feyn- 
man [T5], and later Margolus [15] , suggested that quantising this model was 
important, for two reasons: firstly, because in CA computation occurs without 
extraneous (unnecessary) control, hence eliminating a source of decoherence; 
and secondly because they are a good framework in which to study the quan- 
J^~j turn simulation of a quantum system. From a computation perspective there 

rS are other reasons to study QCA, such as studying space-sensitive problems in 

C3 computer science, e.g. 'machine self-reproduction' [SH] or 'Firing Squad Synchro- 

nisation', which QCA allow in the quantum setting. There is also a theoretical 
physics perspective, where CA are used as toy models of quantum space-time 
[17] . The first approach to defining QCA |2ll5l37j was later superseded by a 
more axiomatic approach [8 9 30] together with the more operational approaches 
[11124128129134157) . 

The most well known CA is Conway's 'Game of Life', a two-dimensional CA 
which has been shown to be universal for computation, in the sense that any 
Turing Machine (TM) can be encoded within its initial state and then executed 
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by evolution of the CA. Because TM have long been regarded as the best def- 
inition of 'what an algorithm is' in classical computer science, this result could 
have been perceived as providing a conclusion to the topic of CA universality. 
This was not the case, because CA do more than just running any algorithm. 
They run distributed algorithms in a distributed manner, model phenomena to- 
gether with their spatial structure, and allow the use of the spatial parallelism 
inherent to the model. These features, modelled by CA and not by TM, are all 
interesting, and so the concept of universality must be revisited in this context 
to account for space. This is achieved by returning to the original meaning of 
the word universality |1I10I12) . namely the ability for one instance of a compu- 
tational model to be able to simulate other instances of the same computational 
model. Intrinsic universality formalises the ability of a CA to simulate another in 
a space-preserving manner |20 |26|32] , and was extended to the quantum setting 
in [51515] . 

There are several related results in the CA literature. For example, refs. 
[18122123] provide computation universal Reversible Partitioned CA construc- 
tions, whereas ref. [21] deals with their ability to simulate any CA in the one- 
dimensional case. The problem of minimal intrinsically universal CA was ad- 
dressed, cf. [27], and for Reversible CA (RCA) the issue was tackled by Durand- 
Lose [13114] . The difficulty is in having an n-dimensional RCA simulate all other 
n-dimensional RCA and not, say, the (n — 1) -dimensional RCA, otherwise a 
history- keeping dimension could be used, as by Toffoli [53]. There are also sev- 
eral other QCA related results. Watrous [38] has proved that QCA are universal 
in the sense of QTM. Shepherd, Franz and Werner [3T] defined a class of QCA 
where the scattering unitary £/; changes at each step i (CCQCA). Universality 
in the circuit-sense has already been achieved by Van Dam [33], Cirac and Voll- 
brecht [55], Nagaj and Wocjan [24 and Raussendorf 29 . In the bounded-size 
configurations case, circuit universality coincides with intrinsic universality, as 
noted by Van Dam |34j . QCA intrinsic universality in the one-dimensional case 
is resolved in ref. 0]. Given the crucial role of this in classical CA theory, the 
issue of intrinsic universality in the n-dimensional case began to be addressed 
in refs. |5l6j . where it was shown that a simple subclass of QCA, namely Parti- 
tioned QCA (PQCA), are intrinsically universal. Having shown that PQCA are 
intrinsically universal, it remains to be shown that there exists a n-dimensional 
PQCA capable of simulating all other n-dimensional PQCA for n > 1, which is 
presented here. 

PQCA are QCA of a particular form, where incoming information is scattered 
by a fixed unitary U before being redistributed. Hence the problem of finding 
an intrinsically universal PQCA reduces to finding some universal scattering 
unitary U (this is made formal in section 2.2 see Fig{2]). Clearly the universality 



requirement on U is much more difficult than just quantum circuit universality. 
This is because the simulation of a QCA H has to be done in a parallel, space- 
preserving manner. Moreover we must simulate not only an iteration of H but 
several (H 2 , . . . ), so after every simulation the universal PQCA must be ready 
for a further iteration. 
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From a computer architecture point of view, this problem can be recast in 
terms of finding some fundamental quantum processing unit which is capable of 
simulating any other network of quantum processing units, in a space-preserving 
manner. From a theoretical physics perspective, this amounts to specifying a 
scattering phenomenon that is capable of simulating any other, again in a space- 
preserving manner. 

2 An Intrinsically Universal QCA 

The aim is to find a particular {/-defined PQCA which is capable of intrinsically 
simulating any ^-defined PQCA, for any V. In order to describe such a U- 
defined PQCA in detail, two things must be given: the dimensionality of the cells 
(including the meaning attached to each of the states they may take), and the 
way the scattering unitary U acts upon these cells. First we discuss the general 
scheme used to solve this problem, and then describe the PQCA implementing 
it. The necessary definitions for rt-dimensional QCA are given in refs. 5 6 . 

2.1 Circuit Universality versus Intrinsic Universality in Higher 
Dimensions 

As already discussed, intrinsic universality refers to the ability for one CA to 
simulate any other CA in a way which preserves the spatial structure of the 
simulated CA. Conversely, computation universality refers to the ability of a CA 
to simulate any TM, and hence run any algorithm. Additionally, circuit univer- 
sality is the ability of one CA to simulate any circuit. These are Nand gate 
circuits for classical circuits and CA, and Toffoli gate circuits for reversible 
circuits and CA. Informally, in a quantum setting, circuit universality is the 
ability of a PQCA to simulate any unitary evolution expressed as a combination 
of a universal set of quantum gates, such as the standard gate set: Cnot, R(f ) 
(also known as the | gate), and the Hadamard gate. The relationships be- 
tween these three concepts of CA universality have been noted previously [T2] . 
A computation universal CA is also a circuit universal CA, because circuits are 
finitary computations. Moreover, an intrinsic universal CA is also a computation 
universal CA, because it can simulate any CA, including computation universal 
CA. Hence intrinsic universality implies computation universality, which implies 
circuit universality. 

In one-dimension this is not an equivalence. Intuitively, computation uni- 
versality requires more than circuit universality, namely the ability to loop the 
computation, which is not trivial for CA. Similarly, intrinsic universality requires 
more than computation universality, such as the ability to simulate multiple com- 
municating TM. In the classical setting there are formal results that distinguish 
these ideas [26] . 

In n-dimensions, it is often assumed in the classical CA literature that circuit 
universality implies intrinsic universality, and hence these are all equivalent [26 . 
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Strictly speaking this is not true. Consider a two-dimensional CA which runs one- 
dimensional CA in parallel. If the one-dimensional CA is circuit /computation 
universal, but not computation/intrinsically universal, then this is also true for 
the two-dimensional CA. Similarly, in the PQCA setting, the two-dimensional 
constructions in |28j and |29j are circuit universal but not intrinsically universal. 

However, this remains a useful intuition: Indeed, CA admit a block represen- 
tation, where these blocks are permutations for reversible CA, while for PQCA 
the blocks are unitary matrices. Thus the evolution of any (reversible/quantum) 
CA can be expressed as an infinite (reversible/quantum) circuit of (reversible/ 
quantum) gates repeating across space. If a CA is circuit universal, and if it 
is possible to wire together different circuit components in different regions of 
space, then the CA can simulate the block representation of any CA, and hence 
can simulate any CA in a way which preserves its spatial structure. It is intrinsi- 
cally universal. This is the route followed next in constructing the intrinsically 
universal n-dimensional PQCA. First the construction of the 'wires', which can 
carry information across different regions of space, is considered. Here these 
are signals which can be redirected or delayed using barriers, with each signal 
holding a qubit of information. Secondly, the 'circuit-pieces' are constructed, by 
implementing quantum gates which can be combined. One and two qubit gates 
are implemented as obstacles to, and interactions of, these signals. 

2.2 Flattening a PQCA into Space 

In the classical CA literature it is considered enough to show that the CA imple- 
ments some wires carrying signals, and some universal gates acting upon them, 
to prove that an n-dimensional CA is in fact intrinsically universal. Any CA can 
be encoded into a 'wire and gates' arrangement following the above argument, 
but this has never been made explicit in the literature. This section makes more 
precise how to flatten any PQCA in space, so that it is simulated by a PQCA 
which implements quantum wires and universal quantum gates. Flattening a 
PQCA means that the infinitely repeating, two-layered circuit is arranged in 
space so that at the beginning all the signals carrying qubits find themselves in 
circuit-pieces which implement a scattering unitary of the first layer, and then 
all synchronously exit and travel to circuit-pieces implementing the scattering 
unitary of the second layer, etc. An algorithm for performing this flattening can 
be provided, however the process will not be described in detail, for clarity and 
following the classical literature, which largely ignores this process. 

The flattening process can be expressed in three steps: Firstly, the ^-defined 
PQCA is expanded in space by coding each cell into a hypercube of 2™ cells. 
This allows enough space for the scattering unitary V to be applied on non- 
overlapping hypercubes of cells, illustrated in the two-dimensional case in Fig.[l] 
Secondly, the hypercubes where V is applied must be connected with wires, 
as shown in Fig. [T] (right). Within these hypercubes wiring is required so that 
incoming signals are bunched together to undergo a circuit implementation of V, 
and are then dispatched appropriately, as shown in Fig. [2j (left). This requires 
both time and space expansions, with factors that depend non-trivially (but 



A Simple n-Dimensional Intrinsically Universal QCA 



5 



O • 
O • 



V 

• 


• v 


v° 


•v 



o 




-ji 






o < 

O « 


) 




° 




K 





Fig. 1. Flattening a PQCA into a simulating PQCA. Left: Consider four cells 
(white, light grey, dark grey, black) of a PQCA having scattering unitary V. The 
first layer PQCA applies V to these four cells, then the second layer applies V 
at the four corners. Right: We need to flatten this so that the two-layers become 
non-overlapping. The first layer corresponds to the centre square, and the second 
layer to the four corner squares. At the beginning the signals (white, light grey, 
dark grey, black) coding for the simulated cells are in the centre square. They 
undergo V, and are directed towards the bottom left, top left, bottom right, and 
top right squares respectively, where they undergo V but paired up with some 
other signals, etc. 



uninterestingly) upon the size of the circuit implementation of V and the way 
the wiring and gates work in the simulating PQCA. Next, an encoding of the 
circuit description of the scattering unitary V is implemented in the simulating 
PQCA upon these incoming bunched wires, as shown in Fig. [2] (right). This 
completes the description of the overall scheme according to which a PQCA 
that is capable of implementing wires and gates is also capable of intrinsically 
simulating any PQCA, and hence any QCA. A particular PQCA that supports 
these wires and gates can now be constructed. 



2.3 Barriers and Signals Carrying Qubits 

Classical CA studies often refer to 'signals' without an explicit definition. In this 
context, a signal refers to the state of a cell which may move to a neighbouring 
cell consistently, from one step to another, by the evolution of the CA. Therefore 
a signal would appear as a line in the space-time diagram of the CA. These 
lines need to be implemented as signal redirections. A 2D solution is presented 
here, but this scheme can easily be extended to higher dimensions. Each cell 
has four possible basis states: empty (e), holding a qubit signal (0 or 1), or a 
barrier (■). The scattering unitary U of the universal PQCA acts on 2 x 2 cell 
neighbourhoods. 

Signals encode qubits which can travel diagonally across the 2D space (NE, 
SE, SW, or NW). Barriers do not move, while signals move in the obvious way 
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Fig. 2. Flattening a PQCA into a simulating PQCA (cont'd). Left: Within the 
central square the incoming signals are bunched together so as to undergo a 
circuit which implements V, and are then dispatched towards the four corners. 
This diagram does not make explicit a number of signal delays, which may be 
needed to ensure that they arrive synchronously at the beginning of the circuit 
implementing V. Right: Within the central rectangle, the circuit which imple- 
ments V is itself a combination of smaller circuits for implementing a universal 
set of quantum gates such as Cnot, Hadamard and the R(f ), together with 
delays. These are implemented as explained in sections [23| and [2T] 



if unobstructed, as there is only one choice for any signal in any square of four 
cells. Hence the basic movements of signals are given by the following four rules: 





s 







where s £ {0, 1} denotes a signal, and blank cells are empty. 

The way to interpret the four above rules in terms of the scattering unitary 



U is just case-by-case definition, i.e. they show that U 

Moreover, each rule can be obtained as a rotation of another, hence by stating 
that the [/-defined PQCA is isotropic the first rule above suffices. This convention 
will be used throughout. 

The ability to redirect signals is achieved by 'bouncing' them off walls con- 
structed from two barriers arranged either horizontally or vertically: 



m 



where s again denotes the signal and the shaded cells denote the barriers which 
causes the signal to change direction. If there is only one barrier present in the 
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four cell square being operated on then the signal simply propagates as normal 
and is not deflected: 

Using only these basic rules of signal propagation and signal reflection from 
barrier walls, signal delay (Fig. [3]) and signal swapping (Fig. [4]) tiles can be 
constructed. All of the rules presented so far are permutations of some of the 
base elements of the vector space generated by 

^u;,a;,y,,s£{e,0,l,M} 

therefore U is indeed unitary on the subspace upon which its action has so far 
been described. 
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Fig. 3. The 'identity circuit' tile, an 8 x 14 tile taking 24 time-steps, made by 
repeatedly bouncing the signal from walls to slow its movement through the tile. 
The dotted line gives the signal trajectory, with the arrow showing the exit point 
and direction of signal propagation. The bold lines show the tile boundary. 



2.4 Gates 

To allow a universal set of gates to be implemented by the PQCA, certain 
combinations of signals and barriers can be assigned special importance. The 
Hadamard operation on a single qubit-carrying signal can be implemented by 
interpreting a signal passing through a diagonally oriented wall, analogous to a 
semitransparcnt barrier in physics. This has the action defined by the following 
rule: 
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Fig. 4. The 'swap circuit' tile, a 16 x 14 tile, where both input signals are per- 
muted and exit synchronously after 24 time-steps. As the hrst signal (bottom 
left) is initially delayed, there is no interaction. 



This implements the Hadamard operation, creating a superposition of configu- 
rations with appropriate phases. Using this construction a Hadamard tile can 
be constructed (Fig. [5]) by simply adding a semitransparent barrier to the end 
of the previously defined delay (identity) tile (Fig. [3]). A way of encoding two 
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Fig. 5. The 'Hadamard gate' tile applies the Hadamard operation to the input 
signal. It is a modification of the identity circuit tile, with a diagonal (semitrans- 
parent) barrier added at the end which performs the Hadamard operation. 



qubit gates in this system is to consider that two signals which cross paths in- 
teract with one another. The controlled-R( j) operation can be implemented by 
considering signals that cross each other as interacting only if they are both 1, 
in which case a global phase of is applied. Otherwise the signals continue as 
normal. This behaviour is defined by the following rule: 
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where {0, 1}. This signal interaction which induces a global phase change 

allows the definition of both a two signal controlled-R(^) tile (Fig. |6j) and a 
single signal R(f ) operation tile (Fig. Uj. These rules are simply a permutation 
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Fig. 6. The 'controlled-R(^) gate' tile applies the controlled-R( ?) operation to 
the two input qubits, by causing the signals to interact at the highlighted point 
(grey circle). The qubits are then synchronised so that they exit at the same 
time along their original paths. No swapping takes place. 



and phase change of base elements of the form: 



X 




y 





(and their rotations), therefore U is a unitary operation on the subspace upon 
which its action has so far been described. Wherever U has not yet been defined, 
it is the identity. Hence U is unitary. 



2.5 Circuits: Combining Gates 

A signal is given an 8 x 14 tile (16 x 14 for two signal operations) in which the 
action is encoded. The signals enter each tile at the fifth cell from the left, and 
propagate diagonally NE. Each time step finds the tile shifted one cell to the 
right to match this diagonal movement, giving a diagonal tile. The signal exits 
the tile 14 cells North and East of where it entered. This allows these tiles to 
be composed in parallel and sequentially with the only other requirement being 
that the signal exits at the appropriate point, i.e. the fifth cell along the tile, 
after 24 time-steps. This ensures that all signals are synchronised as in Fig. [2] 
(right), allowing larger circuits to be built from these elementary tiles by simply 
plugging them together. Non-contiguous gates can also be wired together using 
appropriate wall constructions to redirect and delay signals so that they are 
correctly synchronised. 
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Fig. 7. The 'R(f) gate' tile. This tile makes use of a signal, set to |1), which 
loops inside the grid every six time-steps, ensuring that it will interact with 
the signal that enters the tile, and causing it to act as the control qubit to a 
controlled-R(^) operation. It therefore acts as a phase rotation on the input 
qubit, which passes directly through. After 24 time-steps the auxiliary control 
signal has returned to its origin, unchanged, hence the tile can be reused. 

The implemented set of quantum gates, the identity, Hadamard, swap, R( |) 
and controllcd-R(^), gives a universal set. Indeed the standard set of cNot, H, 
R(f ) can be recovered as follows: 

CNot = (I® J ff)(cR( 7 r/4)) 4 (I <g> H) \tp) 

where cR(^) 4 denotes four applications of the controllcd-R(^) gate, giving the 
controlled-PHASE operation. 

3 Conclusion 

This paper presents a simple PQCA which is capable of simulating all other 
PQCA, preserving the topology of the simulated PQCA. This means that the 
initial configuration and the forward evolution of any PQCA can be encoded 
within the initial configuration of this PQCA, with each simulated cell encoded as 
a group of adjacent cells in the PQCA, i.e. intrinsic simulation. The construction 
in section[2]is given in two-dimensions, which can be seen to generalise to n > 1- 
dimcnsions. The main, formal result of this work can therefore be stated as: 

Claim 1 There exists an n- dimensional U -defined PQCA, G, which is an in- 
trinsically universal PQCA. Let H be a n- dimensional V -defined PQCA such 
that V can be expressed as a quantum circuit C made of gates from the set 
Hadamard, Cnot, and R(f )• Then G is able to intrinsically simulate H. 

Any finite-dimensional unitary V can always be approximated by a circuit 
C(V) with an arbitrary small error e = max^) \\V \ip) — C \ip) ||. Assuming in- 
stead that G simulates the C(V)-defined PQCA, for a region of s cells over a 
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period t, the error with respect to the ^-defined PQCA will be bounded by ste. 
This is due to the general statement that errors in quantum circuits increase, 
at most, proportionally with time and space [25) . Combined with the fact that 
PQCA are universal |5l6j . this means that G is intrinsically universal, up to this 
unavoidable approximation. 



3.1 Discussion and Future Work 



QC research has so far focused on applications for more secure and efficient 
computing, with theoretical physics supporting this work in theoretical com- 
puter science. The results of this interdisciplinary exchange led to the assump- 
tions underlying computer science being revisited, with information theory and 
complexity theory, for example, being reconsidered and redeveloped. However, 
information theory also plays a crucial role in the foundations of theoretical 
physics. These developments are also of interest in theoretical physics studies 
where physical aspects such as particles and matter are considered; computer 
science studies tend to consider only abstract mathematical quantities. Univer- 
sality, among the many computer science concepts, is a simplifying methodology 
in this respect. For example, if the problem being studied crucially involves 
some idea of interaction, universality makes it possible cast it in terms of in- 
formation exchanges together with some universal information processing. This 
paper presents an attempt to export universality as a tool for application in 
theoretical physics, a small step towards the goal of finding and understanding a 
universal physical phenomenon, within some simplified mechanics. Similar to the 
importance of the idea of the spatial arrangement of interactions in physics, in- 
trinsic universality has broader applicability than computation universality and 
must be preferred. In short, if only one physical phenomenon is considered, it 
should be an intrinsically universal physical phenomenon, as it could be used to 
simulate all others. 

The PQCA cell dimension of the simple intrinsically universal construction 
given here is four (empty, a qubit (|0) or |1)), or a barrier). In comparison, the 
simplest classical Partitioned CA has cell dimension two [T5]. Hence, although 
the intrinsically universal PQCA presented here is the simplest found, it is not 
minimal. In fact, one can also manage [7] an intrinsically universal PQCA with 
a cell dimension of three, in two different ways. One way is to encode the spin 
degree of freedom (0 and 1) into a spatial degree of freedom, so that now the 
semitransparent barrier either splits or combines signals. The second way is to 
code barriers as pairs of signals as in the Billiard Ball CA model [TH]. These 
constructions may be minimal, but are not as elegant as the one presented here. 
In future work we will show that there is a simple, greater than two-dimensional 
PQCA which is minimal, as it has a cell dimension of two. 
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